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PL  -  long  period  P  wave 

S,  direct  and  long  period  shear  (S)  waves 
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Sv  -  vertically  polarized  shear  wave 

A  -  epicentral  distance ,  km 

T  -  period  of  wave,  s 
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h  -  focal  depth,  km 

t  -  travel  time  from  source  to  detector  for  an  X  -  type  wave,  s 

-  average  velocity  of  P  wave  in  the  earth’s  crust 

\?2  -  average  velocity  of  P  wave  in  the  earth’s  mantle 

Ax  -  amplitude  of  a  wave  of  the  X  type 

x  -  S  wave  decay  time,  s  t  -  t(Ag  =  max)  -  t(Ag  «  1/2  max) 


n  -  number  of  data  points 
a  -  standard  deviation 

ar¥  -  stress  in  the  r  -  ¥  plane 
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U  -  displacement  vector  with  components  (u,v,w)  in  the  (r,f,z)  directions 
<p  -  scalar  potential 

A  -  vector  potential 

a,b,a',  b",  A,B,  -  constants 


DIFFERENCES  BETWEEN  EXPLOSIONS  and  EARTHQUAKES  * 


I.  Introduction 

Along  with  the  rapid  progress  of  socialist  reconstruction  In  China  has 
come  an  Increasing  use  of  explosives  in  normal  daily  industrial  mining.  There¬ 
fore,  the  problem  of  how  to  accurately  distinguish  between  explosions  and  earth¬ 
quakes  has  already  become  an  important  problem  of  study  among  seismic  activities 
in  mining  areas.  Nearby,  small  explosions  are  relatively  easy  to  distinguish, 
because  the  P  wave  and  S  wave  amplitude  rations  and  the  S  wave  periods  are  much 
different  from  earthquakes  on  the  seismic  record.  The  identification  of  far-off 
nuclear  explosions  has  already  become  an  independent  area  of  study  with  a  great 
amount  of  work  being  done  by  others  outside  the  country*.  This  paper  will  pri¬ 
marily  discuss  the  statistical  differences  between  medium  distance  (100-500  km), 
medium  strength  (2.5  <  <  4.5)  explosion  and  earthquake  records,  and  offer 

methods  of  distinguishing  explosions  based  on  these  differences. 

Records  of  15  comparatively  large  industrial  explosions  recorded  since  1969 
from  China's  Xinan,  Shanxi,  and  Beijing  (Peking)  Seismic  Observatories  are  used 
in  this  paper.  Each  station  was  equipped  with  Wiechert  vertical  component  seis¬ 
mograph  recording  the  periods  from  T*  0.05  s  to  1.0  s. 
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II.  Statistical  differences  between  explosion  and  earthquake  records 

Sometimes  the  distinctive  features  of  an  explosion  are  very  clear  on  a 

seismogram,  with  no  need  for  a  careful  examination  to  conclude  that  the  event  Is 

an  explosion.  Sometimes  an  explosion  has  very  few  differences  from  an  earthquake 

record  and  is  almost  impossible  to  discern  as  an  explosion.  This  clearly  shows 

that  the  differences  between  explosion  and  earthquake  records  have  laregly 

statistical  variations.  Based  on  the  results  from  large  amounts  of  explosion 

data,  we  shall  now  develop  in  the  discussion  the  differences  between  two  types. 

1.  The  relationship  between  t^_p  and  tg_p 

n 

Assuming  that  the  earth’s  crustal  thickness,  H,  is  40  km,  then  the 

travel  times  for  P  and  P  waves  are: 

n 


A 

V. 


+ 


i  -  (Vj/V^2 


Because  explosions  are  set  off  at  the  earth's  sruface  (h  »  0),  and  the  depth  for 

most  earthquakes  is  h  *  15-30  km,  the  propagation  paths  of  the  seismic  waves 

will  not  be  the  same,  so  the  relationship  between  t^_p  and  tg_p  will  also  be 

n 

different.  Figure  1  shows  the  result  for  a  number  of  explosions  and  earthquakes. 

From  this  figure  it  can  be  seen  that  whether  the  event  is  an  explosion  or  an 

earthquake,  t^_p  and  tg_g  show  a  linear  recursive  relationship  resembling  a 
n 

binary  normal  distributed  random  variation. 

To  obtain  the  straight  line  equation  and  its  parameters,  we  have  for  the 


explosions : 


where 


a'  -  7.62,  b  -  0.305,  o  -  0.478,  n  -  29. 
For  the  earthquake  straight  line  and  its  parameters  we  have: 


tP-P  a'  +  b  t§-P  * 
n 


a'  *  4.56,  b  *  0.304,  a  “  1.09,  ■  6.71,  n  “  41, 


-  In 

S-P  n  l  S-P 


(1) 


33.4 


i-1 

These  two  straight  line  equations  have  clear  differences:  the  points  for  the 

explosions  lie  comparatively  close  to  the  straight  line,  while  the  data  points 

for  the  earthquakes  have  more  scatter,  and  the  two  lines  are  separated  by  about 

3  seconds  (along  the  y-axis].  Now  using  the  assumption  of  a  normal  distribution 

2 

(about  the  lines)  in  the  investigative  method  »  we  form  n-2  independent  equations 
for  t: 


T 

a 


A 


a  -  a 


a 


+  b  ‘s-P 


-  a'  + 


btS-P 


Using  a,  b  for  explosions  and  a,  b  for  earthquakes,  we  obtain 

T  -  17.0  T.  -  0.038 

a  b 


It  appeats  evident  from  the  distribtuion  of  t  that,  using  the  parameters  a  and  b 
to  examine  the  two  straight  lines  as  a  single  probability  distribution  gives 

P(|t|;  T()  -  0 

K|t|-  Tb)  -  1.0 

From  this  it  can  be  seen  that  the  differences  between  the  two  straight  lines  ap- 

j 

i 

! 
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pear  primarily  in  parameter  a.  Because  the  relationship  between  t-_p  and 

~n 

t-_^  depends  primarily  on  the  crustal  structure  and  the  focal  depth,  h,  and  the 
influence  of  other  factors  is  small,  the  parameter  a  is  a  basis  for  distinguish¬ 
ing  explosions. 

2.  The  relationship  between  and  Mg 

The  relationship  between  the  near-earthquake  magnitude,  M^»  and  the 
surface  wave  magnitude.  Mg,  used  in  China  is 

Mg  -  1.13  M^  -  1.08 

Gu  Dengbao  and  Li  Vongte's  results  give 

Mg  -  1.27  (M^  —  1)  —  0.016  m£ 

These  equations  clearly  show  that  the  local  magnitude  of  earthquakes  is 
larger  than  the  surface  wave  magnitude  M  .  However,  this  is  not  the  case  with 
explosions.  Using  a  May,  1971,  explosion  as  an  example,  the  results  of  Yunnan 
Province1 s  Jaotong  and  Kunming  Observatories  are  (Mg  determined  from  Jishi 
[Basic]  seismograph,  Mj^  determined  from  Wiechert  seismograph): 

Jaotong:  Mg  -  M^  *  4.7  -  3.8  ■  0.9 

Kunming:  Mg  -  M^  *  4.9  -  3.7  -  1.2 

Mg  is  clearly  larger  than  M^  for  explosions,  exactly  the  opposite  of  what  was 
found  for  earthquakes,  so  the  quantity  Mg-M^^  is  one  of  the  differences  between 
the  two. 

There  has  already  been  a  large  amount  of  work  done  outside  of  China  on  the 

use  of  a  similar  relationship  between  the  body  wave  magnitude,  m,  and  the  surface 

1  3 

wave  magnitude,  M  ,  to  detect  underground  nuclear  explosions  9  • 
b 

3.  The  Sy/Sjj  ratio 

If  an  explosion  is  detonated  in  an  infinite  uniform  space,  then  only  longitu¬ 


dinal  waves,  P,  will  be  produced.  But  actual  explosion  records  show  not 


only  P  waves,  but  also  S  waves.  This  is  because  actual  explosions  are  set  off 
in  a  non-uniform  half-space.  Based  on  the  symnetry  of  an  explosion,  however. 


n 

( 


the  Sg  portion  of  the  S  wave  should  be  very  small.  Starting  out  from  this 
consideration,  the  S  waves  for  three  explosion  events  were  examined,  and  the 
ratios  Sg/Sy  were  calculated.  The  results  are  shown  below: 


Station 

Vs' 

Xichang 

0.19 

Mabian 

0.05 

Tangdan 

0.06 

Panzhihua 

0.28 

Qulnshul 

0.75 

Huoxlan 

0.54 

Puxian 

0.03 

These  stations  were  distributed  in  different  directions  around  the  explosion. 
For  each  station,  Sg  is  less  than  Sy,  and  in  most  cases  it  is  much  less.  The 


average  of  the  values  from  the  seven  stations  is  Sg/Sy 


0.27. 

4 

Below  we  give  a  simple  analysis,  following  the  method  of  M.  Bath  .  If  we 
have  a  displacement  U  ■  7  $+  VxA  in  free  space,  then  $  and  A  will  separately 
satisfy  the  wave  equations: 


2 

a  d> 
at2 


2  2 

*  7^  , 


4 .  bV  * . 

at 


Under  normal  conditions,  an  explosion  shows  at  least  cylindrical  symmetry,  so 
that  in  cylindrical  co-ordinates: 


3A  3A 

r  m _ z  m  3£ 

ay  “  ay  “  ay 


He  now  write  the  r,  y,  and  z  components  of  U  as  u,  v,  and  w  below: 


The  above  equation  shows  that  in  the  r-z  plane  only  $  and  A  are  involved  in  the 

components  of  motion;  these  are  equivalent  to  the  P  and  S  waves.  In  the  ¥ 

direction  only  A  and  A  are  non-zero,  corresponding  to  the  Su  wave.  Because  of 
r  z  n 

the  axial  symmetry  we  can  divide  the  and  components.  We  first  discuss  the 
Sjj  wave: 

3Ar  3Az 
V  "  3z  *  3r 

Because  of  the  axial  symmetry  we  have: 


and  since  r  is  the  time  of  the  cylindrical  radius  of  the  source,  v*0,  therefore 
v(r,z)  *  0* 

From  this  calculation  it  can  be  seen  that  in  cases  of  axial  symmetry  no 
waves  will  exist,  but  actual  explosions  cannot  staisfy  the  axial  or  spherical 
symmetry  requirement  because  of  the  influence  of  the  structure  of  the  medium. 

So  waves  will  appear,  but  will  be  much  smaller  than  Sv  waves.  This  is  a 
clear  difference  between  explosions  and  earthquakes. 

4.  The  ratio  of  Ap/Ag 

The  statistical  results  from  most  earthquakes  shows  that  the  P  wave 
amplitude  is  much  smaller  than  the  S  wave  amplitude,  while  the  results  from  many 
explosions  shows  that  the  P  wave  and  S  wave  amplitudes  are  more  similar,  with 
some  records  showing  a  much  larger  P  than  S  wave  amplitude.  This  difference 
probably  arises  because  of  the  two  different  types  of  source  mechanisms.  Lamb^ 
was  the  earliest  to  calculate  the  ground  motion  produced  by  a  concentrated 
perpendicular  force,  and  in  1914  Sharpe^  calculated  a  spherical  explosive  source 
in  an  infinite  medium.  After  this  there  were  many  calculations  concerning 
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explosive  sources.  All  of  the  calculations  show  from  a  theoretical  standpoint 
the  comparatively  large  amplitude  of  the  P  waves.  From  the  figures  [2  and  3] 


one  can  see  that  the  A  /Ac  ratio  for  earthquakes  has  small  scatter,  with  the 

P  » 

values  falling  between  0  and  1.  For  explosions,  the  scatter  in  the  values  of 
the  ratio  is  larger,  with  the  largest  value  being  2.8.  There  are  many  factors 
influencing  the  A  /Ac  ratio  for  explosions,  the  most  important  factors  include 

p  b 

the  explosive  mechanism,  the  structure  of  the  medium,  the  structure  of  the 

observatory  foundation,  [seismograph  pier,]  etc.  Because  of  these  factors,  a 

single  explosion  record  at  an  observatory  may  have  a  small  A  /A  ratio,  but  from 

p  b 

statistical  observations  at  many  stations  the  ratio  for  explosions  will  be  much 
larger  than  for  earthquakes. 

5.  Average  period,  T,  of  the  S  waves 

In  general  it  can  be  said  that  for  observations  at  the  same  distances, 
the  S  waves  from  explosions  have  a  larger  period  than  those  from  earthquakes. 

The  period  of  S  waves  from  small,  nearby  explosion  records  are  about  0. 1-0.2 
seconds  longer  than  those  from  earthquake  records  (see  figures  4,  5,  and  6). 
Sometimes  on  the  record  of  a  small,  nearby  explosion,  the  P  and  S  waves  both 
exhibit  a  stable  period  of  vibration  [i.e.,  a  steady  monochromatic  frequency 
content]  (figures  7  and  8)  which  is  never  observed  for  earthquakes. 

Explosion  records  in  the  range  of  100-500  km  have  other  distinctive  fea¬ 
tures.  Earthquake  records  have  large  changes  in  period  from  the  P  wave  to  the  S 
wave  tail,  while  the  change  in  the  period  on  explosion  records  is  comparatively 
small.  The  explosion  records  are  smooth  and  sleek  in  appearance,  with  smaller 
high  frequency  perturbations.  Figure  9  is  a  record  of  the  1969  Pohai  earthquake 
recorded  at  Madaoyu  station,  and  figure  10  is  a  record  of  a  1970  explosion  at 
the  same  station.  In  these  two  figures  it  is  very  easy  to  see  these  distinctive 
features. 
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According  to  the  results  of  Nersesov  ,  the  primary  period  of  the  explosion 
depends  mainly  on  the  charge  size,  and  the  distance  from  the  shot  is  not  as 
important.  We  have  examined  the  mean  S  wave  period,  T,  for  a  number  of  ex¬ 
plosion  and  earthquake  records.  Figure  11  shows  the  relationship  between  T  and 
M^.  From  this  figure  it  can  be  seen  that  the  period  of  the  S  waves  from  explosions 
is  0.15  to  0.20  seconds  larger  than  the  period  for  earthquakes. 

The  basic  reason  for  the  longer  explosion  wave  periods  is  the  different 
mechanism  in  explosions  and  earthquakes.  For  explosions,  the  expanding  force  is 

applied  for  a  very  short  period  of  time.  According  to  the  experimental  and  ob- 

9  10 

servational  results  of  Rulev  *  ,  earthquake  records  are  composed  of  P,  P  ,  and 

l j 

S^,  and  waves,  with  the  P^  and  waves  arriving  separately  from  and  after 

the  P  and  S  waves  in  a  long  period  wave  train  with  characteristics  of  both  body 

and  surface  waves.  In  explosion  records,  however,  only  P,  P  ,  and  S  waves  are 

Li  Li 

present,  with  none  of  the  comparatively  short-period  S  waves  (see  figures  9  and 
10),  so  that  the  S  component  waves  on  explosive  records  have  an  average  period 
longer  than  the  S  component  waves  on  earthquake  records. 

6.  S  wave  decay  time  t 

The  expanding  force  during  an  explosion  is  in  effect  for  only  a  very 
short  period  of  time:  on  the  order  of  milliseconds  for  most  nuclear  explosions, 
and  slightly  longer  for  conventional  (industrial)  explosions.  For  earthquake 
sources,  however,  at  least  for  earthquakes  along  linear  faults,  many  people  con¬ 
sider  them  to  be  caused  by  the  process  of  sliding  along  a  fault;  accordingly, 
the  spreading  speed  of  the  break  can  be  calculated^,  and  it  is  found  that  small 
earthquakes  take  a  few  seconds.  The  larger  the  earthquake,  the  longer  the  dura¬ 
tion  of  the  source.  For  example,  the  1960  Zhili  earthquake  lasted  over  200 
seconds.  The  source  duration  is  therefore  shorter  for  explosions  than  for 
earthquakes,  and  the  larger  their  magnitudes,  the  larger  the  difference.  Because 
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of  this,  the  duration  of  explosion  records  is  shorter  than  the  duration  of 
earthquake  records.  In  particular,  since  shear  stress  is  the  major  source  of 
most  earthquakes,  the  difference  in  the  duration  of  the  S  waves  is  very  evident. 
We  define  the  time  for  the  S  wave  amplitude  to  decline  from  its  maximum  value  to 
one-half  the  maximum  value  as  the  S  wave  decay  time  t.  Figure  12  shows  the 
relationship  between  x  and  for  explosions  and  earthquakes.  From  this  figure 
we  can  see  that  t  is  smaller  for  explosions  than  for  earthquakes,  and  that  the 
discrepancy  is  greater  for  larger  magnitudes. 

7.  Polarity  distribution  of  primary  ground  motion 

The  polarity  distribution  of  the  first  motion  of  most  structured 
earthquakes  has  a  comparative  regularity  in  all  four  quadrants.  This  is  not  the 
case  for  explosions,  since  the  theoretical  first  motions  for  explosions  should 
all  be  positive.  Actually,  owing  to  the  complex  structure  of  the  medium,  the 
conditions  at  each  detector  will  be  different,  so  that  explosions  will  also  show 
some  negative  polarity  first  motions,  but  the  distribution  is  not  regular. 

Figure  13  shows  the  polarity  distribution  of  the  first  motion  for  a  1971  explo¬ 
sion.  This  figure  clearly  shows  that  negative  primary  motions  from  an  explosion 
are  a  normal  phenomena 

8.  Surface  waves 

Love  waves  are  S„  component  surface  waves,  and  Rayleigh  waves  are  S,r 
rl  V 

component  surface  waves.  According  to  earlier  theories,  explosions  should  not 
produce  any  Love  waves^.  After  a  =  4.6  explosion  at  a  certain  place  in  May, 
1971,  a  =  5.0  earthquake  occurred  in  the  same  area  in  September,  offering  a 
very  good  opportunity  for  comparison,  as  shown  in  figures  14  and  15.  A  com¬ 
parison  of  the  two  figures  shows: 


Explosion:  There  is  almost  no  Love  wave  at  all,  the  vertical 

component  Rayleigh  wave  becomes  very  large,  and  the 
record  is  generally  smooth* 

Earthquake:  The  level  of  the  Love  wave  grows  very  large,  and  there 
is  a  lot  of  high  frequency  disturbances  on  the  record* 

/  To  sum  up  the  above  discussion,  based  on  dealing  with  a  large  number  of 

explosion  events,  and  summarizing  the  statistical  differences  between  explosion 

and  earthquake  seismograms  mentioned  above,  we  have  proposed  a  number  of  useful 

methods  for  distinguishing  explosions.  In  actual  practice,  one  must  look  at  all 

of  the  conditions,  and  use  a  number  of  techniques  to  be  successful. 

III.  Discussion 

Because  this  paper  is  based  on  a  review  of  seismograms  from  the  Wiechert 
seismograph,  comparisons  [of  the  results  of  this  paper]  with  other  types  of  in¬ 
struments,  especially  long-period  seismographs,  will  only  have  value  as  a 
[general]  reference.  Owing  to  the  bandwidth  limitations  of  the  instruments,  a 
study  of  the  spectral  differences  was  not  conducted. 
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FIGURE  CAPTIONS 

1  -  Relationship  between  t-_p  and  t-_- 

n 

x  -  axis:  tg_-  (seconds)  y  -  axis:  t^p  (seconds) 

n 

upper  line  -  earthquakes  lower  line  -  explosions 

2  -  Relationship  between  Ap/Ag  and  A  for  earthquakes 

x  -  axis  -  A  (km)  y  -  axis:  Ap/Ag 

3  *  same  as  for  2,  except  for  explosions 

4  -  Dahe  station  earthquake  seismogram 

CS-P  ■  3'3-  «L  -  *-6 

5  -  Dahe  station  explosion  seismogram 

CS-P  ■  2'2-  \  •  '-2 

6  -  Dahe  station  explosion  seismogram 

2-3-  »L  '  °-7 

7  -  Dahe  station  explosion  seismogram 

CS-P  ■  2-2>  «L  ■  °-9 

8  -  Panzhihua  station  explosion  seismogram 

‘s-P  "  °-7’  “l  -  1*5 

9  -  Seismogram  of  Pohai  earthquake  of  August  21,  1968 

=  3.4,  t-_-  ”  32.0,  Madayu  station  record 

10  -  Explosion  during  1970 

=  3.2,  tg_p  *  31.0,  Madayu  b ration  record 

11  -  T,  the  average  period  of  the  S  waves 

x  -  axis:  y  -  axis:  T  (seconds) 

upper  line  -  explosions  lower  line  -  earthquakes 

12  -  S  wave  decay  time  x 

x  -  axis:  y  -  axis:  t  (seconds) 

upper  line  -  earthquakes  lower  line  -  explosions 


FIGURE  CAPTIONS,  Cont 
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14  - 


Polarity  distribution  of  first  motions  for  an  explosion  during  1971 
o  -  seismograph  location  *  -  explosion  location 

+,  -  polarity  of  first  motion 

Chengdu  station  Jlshi  (Basic)  model  seismograph  record  of  an 
earthquake  in  September  1971 
M^  -  5.0,  A  -  496  km 
ditto  14  except  explosion  in  May,  1971 
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03  5 

«.r-2.2  Alt-  1.2. 


m  7 

»!-,-=  2.2  Mt-0.9. 


Figures  4-7. 
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w  «*£*•*£* 

ts~f  *=0.7  A/;  =?*  1.5. 


®9  19««  fpg  21 


A/z.  -  3-4.  ts-f  =  32-0  0,q,iIiS6ia*. 


ill.M 


